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1 Introduction 



This diploma thesis analyses static, spherically symmetric perfect fluid so- 
lutions to Einstein's field equations with cosmological constant. New kinds 
of global solutions are described. 

By a global solution one means an inextendible spacetime satisfying the 
Einstein equations with cosmological constant with a perfect fluid source. 
The matter either occupies the whole space or has finite extend. In the 
second case a vacuum solution is joined on as an exterior field. 

Global static fluid ball solutions with finite radius at which the pressure 
vanishes are called stellar models. 

Recent cosmological observations give strong indications for the presence 
of a positive cosmological constant with A < 3 x 10~^^m~^. On the other 
hand Anti-de Sitter spacetimes, having negative cosmological constant, are 
important in the low energy limit of superstring theory. 

Therefore it is interesting to analyse solutions to the field equations with 
cosmological constant representing for example relativistic stars. 

Vanishing cosmological constant A = 

The first static, spherically symmetric perfect fluid solution with constant 
density was already found by Schwarzschild in 1918. 

In spherical symmetry Tolman j25j and Oppenheimer and Volkoff |18j 
reduced the field equations to the well known TOV equation. 

The boundary of stellar models is defined to be where the pressure van- 
ishes. At this surface a vacuum solution is joined on as an exterior field. 
In case of vanishing cosmological constant it is the Schwarzschild solution. 
For very simple equations of state Tolman integrated the TOV equation and 
discussed solutions. Although he already included the cosmological constant 
in his calculations he did not analyse them. He stated that the cosmological 
constant is too small to produce effects. 

Buchdahl assumed the existence of a global static solution, to show 
that the total mass of a fluid ball is bounded by its radius. He showed 
the strict inequality M < (4/9) -R, which holds for fluid balls in which the 
density does not increase outwards. It implies that radii of fluid balls are 
always larger than the black-hole event horizon. 

Geometrical properties of constant density solutions were analysed by 
Stephani |^[^. He showed that they can be embedded in a five dimen- 
sional flat space and that they are conformally flat. The cosmological con- 
stant can easily be included in his calculations by redefining some variables. 
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Conformal flatness of constant density solutions is easily shown with the use 
of new Buchdahl variables, which will be introduced in the fourth chapter. 
It will be seen that the field equations imply the vanishing of the squared 
Weyl tensor. 

Moreover he remarked that constant density solutions without cosmo- 
logical constant have the spatial geometry of a 3-sphere. With cosmological 
constant the spatial part of the metric may also be Euclidean or hyperbolic, 
depending on the choices of constant density and cosmological constant. 

As already said, Buchdahl assumed the existence of a global static solu- 
tion to derive the upper bound on the mass for given radius. Rendall and 
Schmidt ^Hl proved the existence of such a global static solution. Baum- 
garte and Rendall later improved the argument with less assumptions on 
the equation of state and the pressure. This was further improved by Mars, 
Martin-Prats and Senovilla |14j . 

There is a conjecture that asymptotically flat static perfect fluid solutions 
are spherically symmetric. If this is true then this means that spherically 
symmetric solutions are the most general static perfect fluid models. 

Non-vanishing cosmological constant A 7^ 

This diploma thesis proves existence of a global solution for cosmological 
constants satisfying A < Anph. ph denotes the boundary density and is given 
by the equation of state. For a degenerate neutron star one may assume the 
boundary density to be the density of iron. In physical units this leads to 
A < (47rG/c^)pFe ~ 7 X 10~^'^m~^. This condition is much weaker than the 
mentioned cosmological upper bound. 

For larger values the existence cannot be proved with the used argu- 
ments. Finiteness of the radius for a given equation of state is discussed, a 
necessary and a sufficient condition are shown. 

Another aim is to derive an analogous Buchdahl inequality that includes 
the cosmological term. For positive cosmological constants the vacuum re- 
gion of spacetime may contain a cosmological event horizon. With the analo- 
gous Buchdahl inequality it will be proved that stellar models have an upper 
bound given by that event horizon. 

Collins [3 stated that for a fixed equation of state and cosmological 
constant the choice of central pressure and therefore central density does 
not uniquely determine the solution. This is disproved. 

Static perfect fluid solutions with cosmological constant were analysed 
by Kriele ^3] and later by Winter |2H| • Both derived the analogous TOV-A 
equation. The flrst one shows uniqueness of the solution for given pressure 
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and density distributions, which already disproved Colhns [Hj. An analogous 
type of Buchdahl inequality is derived but not discussed in the context of 
upper and lower bounds on radii of stellar objects. Winter integrates the 
TOV-A equation from the boundary inwards to the centre, without proving 
the existence of that boundary. This leads to solutions with non-regular 
centres and is therefore not suitable for discussing stellar models. 

Constant density solutions with cosmological constant were first analysed 
by Weyl |27j. In this remarkable paper these solutions are described for 
different values of the cosmological constant. The different possible spatial 
geometries were already pointed out and a possible coordinate singularity 
was mentioned. 

More than 80 years later Stuchhk analysed these solutions again. 
He integrated the TOV-A equation for possible values of the cosmological 
constant up to the limit A < Airpo, where po denotes this constant density. 
In these cases constant density solutions describe stellar models. The third 
chapter shows that a coordinate singularity occurs if A > 47rpo, already 
mentioned by Weyl fI7\ . 

If the cosmological constant equals this upper bound the pressure van- 
ishes at the mentioned coordinate singularity. In this case it is not possible 
to join the Schwarzschild-de Sitter solution for the vacuum. One has to 
use the Nariai solution ^1 E| to get the metric at the boundary. For 
larger cosmological constant the pressure will vanish after the coordinate 
singularity. The volume of group orbits is decreasing and there one has 
to join the Schwarzschild-de Sitter solution containing the r = singular- 
ity. Increasing the cosmological constant further leads to generalisations of 
the Einstein static universe. These solutions have two regular centres with 
monotonically decreasing or increasing pressure from the first to the second 
centre. Certainly the Einstein cosmos itself is a solution. Another new kind 
describes solutions with a regular centre and increasing divergent pressure. 
In this case the spacetime has a geometrical singularity. These solutions are 
unphysical and therefore not of great interest. 

With the arguments used in this diploma thesis the existence of a global 
solution with given equation of state can only be proved for cosmological 
constants satisfying the bound mentioned above, A < 47r/jft. 

Unfortunately, with this restriction solutions of the new kinds do not 
occur. Thus one cannot prove existence of the new kinds of solutions for a 
prescribed equation of state. 

In addition to what was said before this paper shortly investigates New- 
tonian limits of the Schwarzschild-de Sitter solution and of the TOV-A equa- 
tion. In the first case both possible horizons shrink to a point in the limit. 
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In the second, as in the case without cosmological term, this leads to the 
Euler equation for hydrostatic equiUbrium. Nonetheless both limits contain 
a A-corrected gravitational potential. 
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2 Cosmological TOV equation 



This chapter deals with the Einstein field equations with cosmological con- 
stant in the spherically symmetric and static case. For a generalised Birkhoff 
theorem see (23ll2n]. 

First a perfect fluid is assumed to be the matter source. This directly 
leads to a A-extended Tolman-Oppenheimer-Volkoff equation which will be 
called TOV- A equation. The TOV- A equation together with the mean den- 
sity equation form a system of differential equations. It can easily be inte- 
grated if a constant density is assumed. 

Next a vacuum solution, namely the Schwarzschild-de Sitter solution, 
is derived. It is the unique static, spherically symmetric vacuum solution 
to the field equations with cosmological constant with group orbits having 
non-constant volume. But there is one other static, spherically symmetric 
vacuum solution, the Nariai solution EI- The group orbits of this 
solution have constant volume. The Nariai solution is needed to join interior 
and exterior solution in a special case but will not be derived explicitly. 

Finally in this introductory chapter Newtonian limits are derived. Limits 
of the Schwarzschild-de Sitter solution and the TOV-A equation are shown. 
The second leads to the fundamental equation of Newtonian astrophysics 
with cosmological constant. 

2.1 The Cosmological Constant 

Originally the Einstein field equations read 

= A^V, (2.1) 

where k is the coupling constant. Later Einstein [Jj introduced the cosmo- 
logical constant mainly to get a static cosmological solution. These field 
equations are 

G'^^ + A5^^ = kV. (2.2) 

The effect of A can be seen as a special type of an energy-momentum tensor. 
It acts as an unusual fluid 

A 

Tjl = — :gf,u = {p^ + —)u^,u^ + P^g^v, 

Ki C 

with P^ = —A/k and equation of state P^ = —p^(?. 
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2.2 Remarks on the Newtonian limit 

The metric for a static, spherically symmetric spacetime can be written 

ds^ = -c^e'^'^'-Ut^ + e^^'Ur'^ + r'^idO'^ + sin^^ d(l?). (2.3) 

Appendix B of shows that in front of the sphere metric is no loss 
of generality for a perfect fluid. This is true for vanishing cosmological 
constant. With cosmological term there is one vacuum solution with group 
orbits of constant volume, the mentioned Nariai solution I17j. 

Metric 1)2. 3|) contains the constant c representing the speed of light. De- 
fine A = roughly speaking the limit A — > corresponds to the New- 
tonian theory. But equations could be non-regular in the limit A = 0. 
Therefore, following jHElISI) the static, spherically symmetric metric has to 
be written 

ds-" = -\e^'^^'Ut' + e^^'^dr^ + r^dO", (2.4) 
A 

with field equations 

G^y + XKg^i, = KT^,y, (2.5) 

where dVl^ = dO'^ + s\v?0 dcfP' . The new term Xi/{r) in (|2.4j) ensures regular- 
ity of the energy- momentum conservation in the limit if matter is present. 
The AA term in the field equations (|2.5|) makes sure that the gravitational 
potential is regular in the limit A = 0. 

2.3 TOV-A equation 

Assume a perfect fluid to be the source of the gravitational field. The 
derivation of the analogous Tolman-Oppenheimer-Volkoff |18l I25j equation 
is shown. Consider metric (|2.4j) 

A 

Units where G = 1 will be used, thus k = SvrA^. A = 1 corresponds to 
Einstein's theory of gravitation in geometrised units. 

The field equations ()2.5() for a perfect fluid can be taken from appendix A, 
(|A.4|) - (|A.7|) with (|A.9|) . These are three independent equations, which im- 
ply energy- momentum conservation (|A.10|) . Thus one may either use three 
independent field equations or one uses two field equations and the energy- 
momentum conservation equation. The aim of this section is to derive a 
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system of differential equations for pressure and mean density, defined later. 
Hence it is more effective to do the second. So consider the first two field 
equations and the conservation equation. This gives 

1 gAi.(r) ^(r- re-'^(^)') - Ae^'^^'^) = 87rp(r)e^'^(^) (2.6) 
Ar^ dr \ J 

^ (l + r\v'[r) - 6"^^")) + AAe"(") = M\^P{r)e<-''^ (2.7) 
-^(AP(r)+/,) = P'(r). (2.8) 

These are three independent ordinary differential equations. But there are 
four unknown functions. Mathematically but not physically one of the four 
functions can be chosen freely. From a physical point of view there are two 
possibilities. Either a matter distribution p = p[r) or an equation of state 
p = p{P) is prescribed. The most physical case is to prescribe an equation 
of state. Equation (|2.6|) can easily be integrated. By putting the constant 
of integration equal to zero because of regularity at the centre one gets 

g-a(r) ^ ^ _ ^g^l r ^2^^^^^^ _ ^ A r g2^^_ ^2.9) 

^ Jo Jo 
Now the definition of 'mass up to r' is used 

m{r) =4tt [ s'^p{s)ds. (2.10) 
Jo 

The second integral of (|2.9|) can be evaluated to finally give 

= 1 - A^ - X^r\ (2.11) 
r 3 ^ ^ 

In addition define the mean density up to r by 

, , m(r) 

w{r) = -j^, (2.12) 

then the above metric component takes the form 

g-a(r) ^ ^ _ X2w{ry - X-r^. (2.13) 

3 

Therefore one may write metric (|2.4|) as 

" a" ^ i-X2w{ry-X^r^^ ^ ' 
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The function i/'(r) can be eliminated from 1)2 .Tp and ()2.8() . The second field 
equation (|2.7|1 yields to 

^ ^ ^^^P{r) + w{r) - f ^ 

2 i_A2u;(r)r2-Afr2' ^ ' ' 

This with (j2IHl) implies the TOV-A equation [211 EHl 

p,, , (A47rP(r) + i/;(r)-f) (AP(r)+p(r)) 

-T r = —r 1 . 2.16 

^ ' 1 - \2w{r)r^ - Af r2 ^ ' 

Putting A = in ()2.16() one finds the Tolman-Oppenheimer-Volkoff equa- 
tion, without cosmological term, short TOV equation, see jl8| I25j. 

If an equation of state p = p{P) is given, the conservation equation 
(|2.8j) can be integrated to give 

rP{->-) 2rlP 

jP^Fr ''''' 

where Pc denotes the central pressure. Using the definition of m{r) then 
1)2.12(1 and (|2.16() form an integro-differential system for p(r) and P{r). Dif- 
ferentiating ()2.12|) with respect to r implies 

w'{r) = - (47rp(P(r)) - 3w{r)) . (2.18) 
r 

Therefore given p = p{P) equations ()2.16|) and ()2.18|) are forming a system 
of differential equations in P{r) and w{r). The solution for w{r) together 
with ()2.1H|) gives the function a(r), whereas the solution of P{r) in ()2.17|) 
defines ^{r). 

There is a alternative way of finding the TOV-A equation if the TOV 
equation is given. Looking back at the Einstein field equations with cosmo- 
logical constant (|2.6|) . p.7j) it is found that they may be derived by a simple 
substitution. Putting 

A 

87r 



A8^' 



in the field equations without cosmological term (|2.1() . gives the A term 
properly. Therefore the effective energy-momentum tensor is defined by 

T'^ ■■= T^u - -g^u. (2.19) 
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In addition it should be remarked how iDefj is given, 

A 

WeS = W + -. 

b 

The possibihty of writing the equations with effective values is of interest 
later. 

2.4 Schwarzschild-anti-de Sitter and Schwarzschild-de Sitter 
solution 

A vacuum solution (T^u = 0) to the Einstein field equations with cos- 
mological constant is derived. For a vanishing cosmological constant the 
Schwarzschild solution follows, for vanishing mass the metric gives the de 
Sitter cosmology. Again one uses metric (|'2.4() 

A 

One has to find a solution to the system of equations 

Gf,^ + XAOf^u = 0. 

In the static, spherically symmetric case this system reduces to three equa- 
tions. But there are only two unknown functions. The third equation is 
not independent. Thus it suffices to consider the first two field equations 
(|A.4|) . ()A.5|) . given in appendix A. 

JLe^Kr) r^-a{r)\ _ j^^Xu{r) ^ q (2.20) 

Ar^ ar \ ) 
^ (l + r\v'(r) - e'^W) + AAe'^^'^) = 0. (2.21) 



r 

The first equation gives 

d_ 

dr 

and can be integrated to 



(^r - re-"(^')) - AAr^ = 



r 3 

where 2M is a constant of integration. The additional A in front of the 
constant of integration is to get interior ()2.14|) and exterior metric continuous 
at a boundary. 
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Again using equation (|2.2U|) . calculating the derivative and multiplying 
with —e""^^'^ one finds 

- (e"('^) + ra'(r) - l) + XAr^e"^''^ = 0. 

The second field equation (|2.2ip can be put in the following form 

(^e'^('^) - rXu'ir) - l) - XAr^e"^''^ = 0. 

Adding up both equations gives a'{r) + Ai^'(r) = 0. This can now be ex- 
pressed with logarithms to give 

Ainfe'^MeMO) =0. 
ar \ / 

One integrates this equation and may set the constant of integration equal 
to one because the time coordinate can be rescaled. This leads to 

and implies the static, spherically symmetric vacuum solution to the field 
equations with cosmological constant, see It is the unique solution 

which group orbits do not have constant volume. 
With the above the following metric is obtained 

d.^ = -1 (l - - A|r^) dt^ + + r'dn\ (2.22) 

Metric ()2.22() is usually called Schwarzschild-de Sitter metric, although it 
was first published by Kottler (1918). It is well defined for radii satisfying 
9''{r) > 0. 

Figure EHl shows Penrose-Carter diagram for the Schwarzschild-de Sitter 
space with A = 1. The constants M and A are chosen such that 9AM^ < 1, 
which means that there are two horizons, see |in( If 9AAf^ = 1 then 
there is only one horizon and if 9AM^ > 1 the spacetime does not contain 
a horizon. 

Figure shows Penrose-Carter diagram for Schwarzschild-anti-de Sitter 
space with A = 1. For positive mass there is always the black- hole event 
horizon. 

Setting A = in ()2.22|) gives the usual Schwarzschild metric. For a 
vanishing mass the de Sitter metric is obtained. 
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2.5 Newtonian limits 



The Newtonian limits of the Schwarzschild-anti-de Sitter and Schwarzschild- 
de Sitter metric and of the TOV-A equation are derived. The frame theory 
approach by Ehlers [U El El will be used to define a one-parameter family of 
general relativistic spacetime models with Newtonian limit. Including the 
cosmological term gives a Newtonian limit with a A-corrected gravitational 
potential. 



2.5.1 Limit of Schwarzschild-anti-de Sitter and Schwarzschild-de 
Sitter models 



A parametrisation of the Schwarzschild-de Sitter metric is given by (|2.22)1 



ds^ = --{1- A A-r^ dt^ + ^ j— + r'dn\ 

In the limit A = the family (A) = |-A5a/3(A),5°^(A),r^^(A)| converges 
to a field of a mass point at the origin. The black-hole event horizon and the 
cosmological event horizon both shrink to a point. The gravitational field is 
given by 



r[JA = 0] = ^-^r, (2.23) 
3 



all others vanish. This corresponds to Newton's equation with cosmological 
term 



N M A 



2.5.2 Limit of the TOV-A equation 

The parametrised metric ()2.14|) is given by 

^ _ 1 A.(r-)^^2 ^2^^2_ 

A 1 - X2w{r)r^ - X^r^ 

The parametrisation of the TOV-A equation 1)2. 16() is 

(4^AP(r) + w{r) - f) (AP(r) + p(r)) 



P'(r) 



1 - 2Xw{r)r^ - X^r^ 
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This differential equation depends on A, it is regular in the limit A = 0. 
Solutions lead to a family of models with Newtonian limit. 

The parametrised TOV-A equation leads to the "fundamental equation 
of Newtonian astrophysics" [22] with cosmological term 



-r^P'ir)[X = 0] = p(r) ^m(r) - 
The gravitational field reads 

r[JA = 0] = ^[A = 0], (2.24) 

all others vanish. Using (|2.15|) in (|2.24|) reproduces (|2.23jl . But with the 
energy-momentum conservation 1)2. 8() it gives 

P'(r) 

r[JA = 0] = -^, (2.25) 

which is the Euler equation for hydrostatic equilibrium in spherical symme- 
try, usually written as 

p{r)f = VP{r). 
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3 Solutions with constant density 



For practical reasons the notation is changed to geometrised units where 
c2 = l/A = 1. 

Assume a positive constant density distribution p = po = const. Then 
w = ^po gives (|2.16() in the form 




In the following all solutions to the above differential equation are de- 
rived, see 1^1231^. In the equation was integrated but not dis- 
cussed for the different possible values of A. In all cases with A < 47rpo 
were discussed using dimensionless variables. This complicated the physi- 
cal interpretation of the results and was an unnecessary restriction to the 
cosmological constant. 

The central pressure Pc = P{r = 0) is always assumed to be positive 
and finite. Using that the density is constant in (|2.17j) metric (|2.14|) reads 

ds^ = -( ^^y^ ) ' dt' + + r^dn^. (3.2) 

\P{r)+poJ i_(i^^o + A)^2 ^ > 

The metric is well defined for radial coordinates r G [0, r) if A > —SirpQ, 
where r denotes the zero of g^^. If A < — Stt/Oq the metric is well defined for 
all r. 

Solutions of differential equation (|3.1|) are uniquely determined by the 
three parameters pQ, Pc and A. Therefore one has a 3-parameter family of 
solutions. 

The two cases of negative and positive cosmological constant are dis- 
cussed separately. In the first case there is no cosmological event horizon 
whereas in the second there is. Analogous Buchdahl |2] inequalities ()3.38|) 
will be derived for the different cases. 

Figures of the pressure are shown. The constant density po a-nd the 
central pressure Pc are equal in all figures, po = Pc = l- Only the cosmo- 
logical constant varies. This follows the approach of the chapter and helps 
to compare the different functions. 

The end of the chapter summarises the different solutions. An overview 
is given. 
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3.1 Spatial geometry of solutions 

Before starting to solve the TOV-A equation one should have a closer look 
at metric ()3.2() . 

The aim of this section is to describe the spatial geometry of metric 
H3.2|) because it depends on the choice of the constant density po and the 
cosmological constant A. For vanishing cosmological constant this was done 
in [Hj. 

For spherically symmetric spacetimes there exists an invariantly defined 
mass function the quasilocal mass. 

For metric (|.S.2|1 quasilocal mass |22] is given by 

, , 8tt o a o , . 

2?nq(r) = — por' + -rl (3.3) 

If R denotes the boundary of a stellar object then 

2mq(i?) = 2M + ^R^, 

where M denotes the total mass of the object. Let 

Svr A 

then the the spatial part of metric (|3.2j) reads 

da^ = ^+rW, (3.4) 

and the quasilocal mass (|3.3() is given by 

2?nq(r) = kr^, (3.5) 

For all positive values of k spatial metric (|3.4)) describes one half of a 
3-sphere of radius 1/^/k. The quasilocal mass mq(r) is positive. Introducing 
a new coordinate by 

r = —= sin a, (3.6) 



gives 



da'^ = y{da^ + sin^a dfi^). 
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This shows that ()3.2() has a coordinate singularity. The metric with radial 
coordinate r only describes one half of the 3-sphere. Geometrically the coor- 
dinate singularity is easily explained. The volume of group orbits r = const, 
has an extremum at the equator r = 1/ y/k. Thus this coordinate cannot be 
used to describe the second half of the 3-sphere. 

Differential equation 1)3.11) is singular at f = 1/^/k. With the new 
coordinate a the differential equation is regular at the the corresponding 
a(r) = 7r/2. 

If A; = then mq(r) = and the spatial metric is purely Euclidean in 
spherical coordinates. In Cartesian coordinates this simply reads 

da'^ = dx^ + dy"^ + dz"^. 

Finally, negative values of k in the spacial metric correspond to three- 
dimensional hyperboloids with negative quasilocal mass. Then using 

r=--^=sinha, (3-7) 

gives 

da'^ = - — —{da^ + sinh^a dQ"^). 
{-k) 

The coordinate singularity becomes important if A exceeds an upper 
limit. The restriction to A < 47rpo in [2l| is due to this fact. Up to this 
limit the radial coordinate r and the new coordinate a are used side by side. 
Although the radial coordinate is less elegant the physical picture is clearer. 

In appendix B it will be shown that constant density solutions are con- 
formally flat |2j . Conformally flat constant density solutions were analysed 
in [2ni and were called generalised solutions of the interior Schwarzschild 
solution. 

If the density is decreasing or increasing outwards solutions are not con- 
formally flat. 

It will also be seen that models with a real geometric singularity exists. 
They are physically not relevant because the pressure is divergent at the 
singularity. 

3.2 Solutions with negative cosmological constant 

In the Schwarzschild-anti-de Sitter spacetime there exists a black-hole event 
horizon. It will be shown that radii of stellar objects are always larger than 
this radius corresponding to the black-hole event horizon. 
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3.2.1 Stellar models with spatially hyperbolic geometry 

A < -Svrpo 

If A < — Svrpo then k and mq(r) are negative. With (|l-i.7j) the denominator 
of can be written as (cosh^a) and the differential equation does not 
have a singularity. 

The volume of group orbits has no extrema, thus metric ()3.2() has no 
coordinate singularities and is well defined for all r. 

Write the differential equation (|3.1() as 



d 



In 



3P + po - A 
P + Po 



Integration gives 



3P(a) + po 



A_ 



P{a) + Po 



d [In cosh a] . 



C cosh a. 



C is the constant of integration, evaluated by defining P{a = 0) 
the central pressure. At the centre it is found that 



C 



3P, + PO - 

Pc + PO 



(3.8) 

(3.9) 

Pc to be 

(3.10) 



Throughout this calculation it is assumed that 87rpo < —A. With Pc > 
this implies that C > 3 and (1 — A/47rpo) > 3, which can easily be verified. 
The constant C will be used instead of its explicit expression given by 1)3. 1U() . 
Then the pressure reads 



1 



P(q) = PO- 



Awpo 



C cosh a 



C cosh a — 3 



(3.11) 



The function P{a) is well defined and monotonically decreasing for all a. 

The pressure converges to — po as a or as the radius r tends to infinity. 
Thus a radius R, where the pressure vanishes, always exists. It is given by 



|87rpo + A| C72 



4ii 



A 



47rpo 



1 



Writing it in q, where ab corresponds to R 



1 



cosh Ob = — 



A 



4vrpo 



(3.12) 



(3.13) 
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shows that the new coordinate a simphfies the expressions . Therefore the 
radial coordinate will only be used to give a physical picture or if it simplifies 
the calculations. 

Prom equation (|3.13|) one can deduce the inverse function. Then one has 
the central pressure as a function of at, Pc = -Pc(«b)- To do this the explicit 
expression for C from 1)3.10(1 is needed. One finds 



1-4^) (coshafe-1) 

Pc = Po^. '-^^ . (3.14) 

(l-i^F^j -3coshab 

The central pressure given by (|3.14() should be finite. Therefore one 
obtains an analogue of the Buchdahl inequality. It reads 

1 A \ 

4vrpo 

Thus there exists an upper bound for aj, for given A and po- Use that 



cosh ah < - ( 1 — — — J • (3.15) 



sinh(arccosh(a)) = \/a'^ — 1, 
then the corresponding radius R is given by 



1 



< (3.16) 

Since the cosmological constant is negative the right-hand side of (|3.16|) is 
well defined. Using the definition of mass M = (47r/3)po-R^ one can rewrite 
(|3.16|) in terms of M, R and A. This leads to 



Uk'^R + R^^-^R^. (3.17) 

This is the wanted analogue of the Buchdahl inequality ()3.38(1 . 

The explicit derivation of the last equation is shown in section ElHl where 
the generalised Buchdahl inequality is derived. 

To emphasise that the cosmological constant is negative write ()3.17|) 



At r = R, where the pressure vanishes, the Schwarzschild-anti-de Sitter 
solution (|2.22|) is joined. Joining interior and exterior solution is needed in 
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Figure 1: Pressure function P{a) with pQ = 1, = 1 and A = — lOvr 

object's radius 




Figure 2: Radius as a function of Figure 3: Radial coordinate 
central pressure R = R{Pc) r = r{a) 

this and many of the following cases. Therefore the next section shows this 
procedure explicitly. 

Figure n shows the typical behaviour of function P{a) Figures 121 

and 131 show the radius of the object as a function of the central pressure 
(|3.12j) and the radius as a function of the new variable a (|3.7j) , respectively. 
Constant density and central pressure are both set to one, po = 1, Pc = 1. 

3.2.2 Joining interior and exterior solution 

At the P = surface the Schwarzschild-anti-de Sitter solution (|2.22|) is 
joined. Since this is needed for all solutions describing stellar objects the 
procedure of joining interior and exterior solution is shown explicitly. 

Gauss coordinates to the r = const, hypersurfaces in the interior are 
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given by X = / Then the interior metric l!\3.'2\\ reads 



2_ / -Pc + Po \ ,^2 , . 2 , /sinhV^x\ ,o2 



The exterior metric, or the vacuum part of the solution, is given by H2.22() . 
Gauss coordinates to the r = const, hypersurfaces are defined by 

dr 

dx 



2A/ _ A 2 
r 3 ' 

The functions x = x(^) ^ = '^(x) s-i'g not given explicitly because they 
involve elliptic functions. However, one does not need the explicit form to 
show that the metric is at Xfei where the pressure vanishes. In Gauss 
coordinates the exterior metric becomes 

= - (i - ^ - yix?) df + dx' + rUfdn'. (3.18) 

Both metrics are joined at P(x = Xb) = 0. Continuity of the metric at Xb 
implies that 

sinh^/^Xfe /X /oim 
= r(xfe). (3.19) 



To get the ^tt-component continuous one has to rescale the time in the 
interior. With 



PO 



both metrics are joined continuously at Xb- Continuity of the first derivatives 
imply 



dr . . / 2M A , , , , 

coshV^Xb = ^iXb) = f-^^- r^Xb). (3.20) 

If one shows that the derivatives of the gt^-components are continuous at Xb 
one is done. The interior part implies 

-iXb) = 21 — - - -r(xfc) 



dx \ r{xb) 3 J Po 
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where P'{xb) can be derived from (|3.1() to give 



dP 



1 sinh V—kxb ( 47r 



yj—k cosh y—kxb 
With (|:-{.2()j) both equations combine to 



A 



do's" 



iXb 



sinh \J —kxb 



dx V-k 
For the exterior metric one finds 



cosh - 



-kxi 



dg^ 

dx 



iXb 



Mxb 



M 



riXbY' 



(A-K A 

I - 3 



2M A 



(3.21) 



r{Xb) 



r{Xb)- 



Use (|3.19|) for the first term and (|3.2nj) for the last term. Furthermore the 
mass is given by 



47r 

M = —po 



sinh - 



-kXb 



"3 



^PoriXbf, 



see (|2.1U() . Then this derivative evaluated at Xb reads 



^sinh V—kXb 



cosh V—kxb 



-Po 



(3.22) 



which equals equation (|3.2H) . Thus the metric is at Xb- 

Since the density is not continuous at the boundary the Ricci tensor is 
not, either. Therefore the metric is at most C^. This cannot be improved. 



3.2.3 Stellar models with spatially Euclidean geometry 

A = -Svrpo 

Assume that cosmological constant and constant density are chosen such 
that Snpo = —A. Then A; = 0, mq{r) = and the denominator of (|3.H) 
becomes one. The differential equation simplifies to 

dP 

— = -ATTr {P + pof , (3.23) 

and the t = const, hypersurfaces of 1)3. 2|) are purely Euclidean. As in the 
former case metric (|3.2I) is well defined for all r. 
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Equation H3.23() can be integrated. If the constant of integration is fixed 
at the centre by P{r = 0) = Pc one obtains 



P{r) 



1 



(3.24) 



27rr2 + 



1 



- Po- 



Pc+PO 



The denominator of the pressure distribution cannot vanish because central 
pressure and density are assumed to be positive. Therefore ()3.24|1 has no 
singularities. 

As the radius tends to infinity the fraction tends to zero and thus the 
pressure converges to —po- 

This implies that there always exits a radius R where P{r = R) = 0. 
Therefore all solutions to ()3.23|1 in ()3.2() describe stellar objects. Their radius 
is given by 



One finds that the radius R is bounded by l/^/Svrpo. Inserting the 
definition of the density yields to 

P 1 1 AttR^ 

R < = 3.26 

27rpo 27r 3M ' ^ ' 

which implies 



Thus (|3.27|) is the analogous Buchdahl inequality to equation (|3.38j) and 
equals (11071) with A = -Svrpo = -6M/R^. 

At r = R, where the pressure vanishes, the Schwarzschild-anti-de Sitter 
solution (|2.22|1 is joined uniquely by the same procedure as before. Since 
the density is not continuous at the boundary the Ricci tensor is not, either. 
Thus the metric is again at most C^. 

Figure 13 shows pressure 1)3. 24() with spatially Euclidean geometry. The 
radius R where the pressure vanishes as a function of central pressure Pc 
(|3.25() is show in figure |SJ As before po = 1, Pc = 1- 




(3.25) 



3M < 2R. 



(3.27) 
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_ central 
5 pressure 



Figure 4: Pressure function P{r) Figure 5: Radius as a function of 
with pq = I, Pc = 1 and A = — Svr central pressure R = R{Pc) 

3.2.4 Stellar models with spatially spherical geometry 

A > -Svrpo 

If A > —Snpo then k and mq{r) are positive. Equation has a singularity 
at 

1 



Stt ^ I A 
—PO + 3 



(3.28) 



In a the singularity f corresponds to a(f) = tt/2. 

As already mentioned in section 13.11 the spatial part of metric (|3.2() now 
describes part of a 3-sphere of radius 1/^/k. The metric is well defined for 
radii less than f. 

Integration of (|3.1|) gives 



-r" 1 + Ccosa 

P(a)=po -'^^ / . (3.29) 



3 — C cos Q 

One finds 



which is less than zero if A < At^pq, which is the restriction in [24] . Since 
this is the considered case the singularity of the pressure gradient is not 
important yet. 

C is a constant of integration and can be expressed by the central pres- 
sure Pc- Again one obtains 

3Pe + Po-i _ (3.31) 

Pc + PO 
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It is the same constant of integration as in ()3.1U() . But possible value for C 
are different. The considered values of the cosmological constant imply that 
1 < C < 3. The pressure is decreasing for all a. 

It was already stated that P['k/2) < 0. Thus there exists an Of, such 
that P{ab) = 0. ab can be derived from equation (|3.29j) and one obtains the 
analogue of (|3.12|1 

cos ab = ^(l - -7—] ■ (3.32) 



C V 47rpo 

It remains to derive the analogue of the Buchdahl inequality for this 
case. One uses (|3.32|) to find Pc = Pc{oib) 



Pc = PO^ T^, (3.33) 

3 cos at — I 1 



47rpo 

which is similar to (|3.12|) . Finiteness of the central pressure in (|3.33|) gives 
With 



sin(arccos(a)) = \/l 
this yields to the wanted analogue 

M4 



< '\ (3.35) 

which equals 1)3. 16(1 . The greater sign reversed because (arccos a) is decreas- 
ing. Thus it can also be rewritten to give 



3M <^R + R^^-^R^, (3.36) 

which equals 1)3. 17(1 . To emphasise that the cosmological constant is negative 
one may write +|A| rather than —A. 

At a = ab 01 r = R the Schwarzschild-anti-de Sitter solution is joined 
and the metric is C^. Without further assumptions this cannot be improved. 

The pressure ()3.29|) is shown in figure Figures [7| and |H1 show the radius 
of the object as a function of the central pressure ()3.32|) and the radius as a 
function of the new variable a (|3.6j) . respectively. 
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Figure 6: Pressure function -P(a) with po = 1, Pc = 1 and A = — 47r 
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Figure 7: Radius clS Si function of Figure 8: 
central pressure R = R{Pc) r = r{a) 



Radial coordinate 
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3.3 Solutions with vanishing cosmological constant 



Assume a vanishing cosmological constant. Then one can use all equations 
of the former case with A = 0. Only one of these relations will be shown, 
namely the Buchdahl inequality PI- It is derived from (|:-{.H5j) 



< (3.37) 
3vrpo 



usmg M = (47r/3)/9o^^ leads to 

4 

M < -R. (3.38) 

Because of the analogous inequalities of the former cases write 

3M < ^R. (3.39) 

For vanishing cosmological constant the pressure is shown in figure IHl Fig- 
ures ^1 and ^2 show the radius of the object as a function of the central 
pressure and the radius as a function of the new variable a (|3.6|1 . respec- 
tively. The plotted pressure and central pressure are given by 1)3. 29() and 
(l3:32ll with A = 0. 




Figure 9: Pressure function P{a) with po = 1, Pc = 1 and A = 
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Figure 10: Radius function of 
central pressure R = R{Pc) 



Figure 11: Radial coordinate 

r = r{a) 



3.4 Solutions with positive cosmological constant 

In the Schwarzschild-de Sitter spacetime there may exist a black-hole event 
horizon and there may also exist a cosmological event horizon. It depends 
on M and A in the Schwarzschild-de Sitter metric 1)2.22(1 which of the cases 
occur. The three possible cases were mentioned in the end of section IT^ 



3.4.1 Stellar models with spatially spherical geometry 

< A < 47rpo 

Integration of the TOV-A equation ()3.1() leads to 



P(a) = po 



47rpo 



1 1 + C COS a 



3 — C cos a 



(3.40) 



Because of (|3.30() the boundary P{R) = exists. As in section 13.2.41 one 
finds 



i (4- ^ 

^2 ^ 3 V 47rpo 



47rpi 







and written in terms of M, i? and A again gives 



2 /4 A 

3M<-R + Rf^--R^. 



(3.41) 



(3.42) 



Since the cosmological constant is positive the square root term is well de- 
fined if 



R < 



4 1 



3 VA 



(3.43) 
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In this section < A < 47rpo is assumed. Using the definition of mass this 
can be rewritten to give 




This leads to 




) 



3 



(3.44) 



where the right-hand side of (|3.44j) has an upper bound given by (|3.42j) . 
If in addition R < l/\/A then the right-hand side of (|3.44|) is bounded by 
one and the exterior Schwarzschild-de Sitter spacetime has two horizons, see 
section [TH Without the additional assumption the exterior spacetime may 
also have one or may also have no horizon. 

The Schwarzschild-de Sitter solution is joined at r = R, where P{R) = 0. 
Defining M = {A7r/3)pQR^ and rescaling the time in the interior the metric 
will be continuous at R. Again Gauss coordinates relative to the hypersur- 
face P{R) = can be used to get the metric at R. Since the density is 
not continuous at the boundary this cannot be improved. 

One can construct solutions without singularities. The group orbits are 
increasing up to R. The boundary of the stellar object r = R can be put 
in region / of Penrose-Carter figure EHl where the time-like killing vector 
is future directed. This leads to figure 1301 the spacetime still contains an 
infinite sequence of singularities r = and space-like infinities r = oo. It 
is possible to put a second object with boundary r = R in region IV of 
Penrose-Carter figure IHDl where the time-like killing vector is past directed. 
This leads to figure |^ and shows that there are no singularities in the 
spacetime. This spacetime is not globally static because of the remaining 
dynamical parts of the the Schwarzschild-de Sitter spacetime, regions He 
and I He- Penrose-Carter diagrams are in appendix C. 

Similarly to the last figures pressure 1)3. 40(1 . central pressure 1)3. 32(1 and 
the new variable a ()3.6)) are shown in figure El figure El and figure ITU 
respectively. 



29 




object' 


radius 


0.3 




0.25 




0.2 




0.15 




0.1 




0.05 





0.5 



1.5 



2.5 



_ central 
3 pressure 




0.2 0.4 0.6 0.8 1 1.2 1.4 



alpha 



Figure 13: Radius function of 
central pressure R = R{Pc) 



Figure 14: Radial coordinate 
r = r(a) 
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3.4.2 Solutions with exterior Nariai metric 



A = Airpo 

In this special case integration of gives 

\ cosa 

^(a = Poii — 7^ • 3.45) 

3 — C cos Q 

The pressure (|3.45|) vanishes at a = = -7r/2, the coordinate singularity in 
the radial coordinate r. 

One would like to join the Schwarzschild-de Sitter metric. But with 
A = 4:7rpo and M = (47r/3)po-R^ it reads 



d.^ = -(l-^)rft^ + ^ + rW. (3.46) 



I _ 3M_ 



The corresponding radius to ab is R = 3M. Therefore one would join both 
solutions at the horizon R = 3M. Since A = 47r/>o the radius R is also given 
by the cosmological constant R = 1/\/A. Thus the black-hole event horizon 
3M and the cosmological event horizon l/\/A are both located at R. 
The interior metric reads 

The volume of group orbits of metric 1)3.46(1 is increasing whereas the group 
orbits of (|3.47|) have constant volume. Therefore it is not possible to join 
the vacuum solution ()3.46|) on as an exterior field to get the metric at 
the P = surface. 

But there is the other spherically symmetric vacuum solution to the 
Einstein field equations with cosmological constant, the Nariai solution jl6l 
I17j . mentioned in chapter El Its metric is 



ds 



2 



1 



-{A cos logr + B sin log rfdt'^ + \dr'^ + dQ'^ 



(3.48) 



A 

where A and B are arbitrary constants. With r = e" this becomes 

ds^ = ^[-{Acosa + B sin a)'^dt'^ + da^ + dn^] . (3.49) 

Metrics (|3.47|) and ()3.49(1 can be joined by fixing the constants A and B. 
With 

a = -^Va, b = ^^±^Va, 



Po Po 
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the metric is at a = 7r/2. As before, since the density is not continuous, 
the metric is at most C^. 

The three figures show P{a), R = R{Pc) and r = r(a), respectively. 
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Figure 15: Pressure function P{a) with pQ = \, Pc = 1 and A = 47r 
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Figure 16: Radius function of 
central pressure R = R{Pc) 



Figure 17: Radial coordinate 
r = r(a) 
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3.4.3 Solutions with decreasing group orbits at the boundary 

47rpo < A < Ao 

Integration of (IH.lj) gives the pressure 



Assume that the pressure vanishes before the second centre a = vr of the 
3-sphere is reached. The condition P{a = vr) < leads to an upper bound 
for the cosmological constant. This bound is given by 



Then 4:Trpo < A < Aq implies the following: 

The pressure is decreasing near the centre and vanishes for some ai,, 
where tt/2 < ai, < tt. Equations (|3.4|) and (|3.6|) imply that the volume of 
group orbits is decreasing if a > 7r/2. 

At Ob one uniquely joins the Schwarzschild-de Sitter solution by M = 
{4:ir/3)poR^ . With Gauss coordinates relative to the P{ab) = hypersurface 
the metric will be . But there is a crucial difference to the former case with 
exterior Schwarzschild-de Sitter solution. Because of the decreasing group 
orbits at the boundary there is still the singularity r = in the vacuum 
spacetime. Penrose-Carter diagram 1321 shows this interesting solution. 

Figure El shows pressure (|3.5U|) . Figures El and EOl show the radius of 
the object as a function of the central pressure and the radius as a function 
of the new variable a, respectively. Constant density and central pressure 
are both set to one, po = 1, Pc = 1. 




(3.50) 




(3.51) 
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Figure 19: Radius function of 
central pressure R = R{Pc) 
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3.4.4 Decreasing solutions with two regular centres 

Ao < A < Ai? 

As before the pressure is given by 



A 



1 ) + C cos a 



3 — C cos Q 

Assume that the pressure is decreasing near the first centre a = 0. This 
gives an upper bound of the cosmological constant 

Aij :=47r/Jo(3Pc/po + l), (3.52) 

where Ag is the cosmological constant of the Einstein static universe. These 
possible values Aq < A < A^; imply: 

The pressure is decreasing is near the first centre a = but remains 
positive for all a because A > Aq. Therefore there exists a second centre at 
a = TT. At the second centre of the 3-sphere the pressure becomes 



1 1 — C 

Pia = vr) = Po ^'"^°3^^ . (3.53) 

It only vanishes if A = Aq. The solution is inextendible. The second centre 
is also regular. This is easily shown with Gauss coordinates. Recall metric 
(|3.2j) written in a 

ds^ = -( \t^ + r (da' + sin2 adi}') . 

\P{a)+poJ ' 

So one already used nearly Gauss coordinates up to a rescaling. With 
X = a/ Vk this becomes 

/ \ 2 

2 f Pc + po Yr,2 , , 2 , |sm^/kx\ ,„2 



Thus the radius in terms of the Gauss coordinate is 

r{x) = -^s,in{Vkx), (3-55) 
and both centres are regular because 

~;~^(x) = cos(Vkx) = =tl ^or Vkx = Oj'?'"- (3.56) 

dx 
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Solutions of this kind are generalisations of the Einstein static universe. 
These 3-spheres have a homogenous density but do not have constant pres- 
sure. They have a given central pressure Pc at the first regular centre which 
decreases monotonically towards the second regular centre. Generalisations 
of the Einstein static universe have been published earlier in Ref. |12j.^ 

The pressure is plotted in figure with po = 1 and Pc = 1 . Figure 
shows the radius function of a. 
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Figure 22: 
r = r(a) 



Radial coordinate 



with /9o = 1, -Pc = 1 and A = 127r 



I would like to thank Aysel Karafistan for bringing this reference to my attention. 
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3.4.5 The Einstein static universe 

A = Ae 

Assume a constant pressure function. Then the pressure gradient vanishes 
and therefore the right-hand side of H3.1|) has to vanish. This gives an unique 
relation between density, central pressure and cosmological constant. One 
obtains (jS^^ 

A = As = 47r {3Pe + po) , 

where Pe = Pc was used to emphasise that the given central pressure cor- 
responds to the Einstein static universe and is homogenous. Its metric is 
given by 

ds"^ = -dt^ + TIe [da^ + sin^ adVt^) , 

where 



V^47r(PE + pe) ' 

TZe 'is the radius of the three dimensional hyper-sphere t = const. The time 
component was rescaled to 1. 

Thus for a given density po there exists for every choice of a central 
pressure Pc a unique cosmological constant given by ()3.52p such that an 
Einstein static universe is the solution to (|3.1j) . 

The constant pressure of the Einstein static universe is plotted in fig- 
ure EHl with pq = 1 and Pc = 1. Figure 121 shows the radius as a function of 
a. Solutions with two centres all have the same radial coordinate r = r(a). 

Pressure radius 
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3.4.6 Increasing solutions with two regular centres 

As < A < A5 

Integration of (|3.1j) again leads to 



= Po ' , / • 3.57 

6 — C COS a 

Assume that the pressure is finite at the second centre. This leads to an 
upper bound of the cosmological constant defined by 

A5 :=47rpo(6Pc/po + 4). (3.58) 

The possible values of the cosmological constant imply: 

The pressure P{a) is increasing near the first regular centre. It increases 
monotonically up to a = tt, where one has a second regular centre. This 
situation is similar to the case where Aq < A < A^;. These solutions are also 
describing generalisations of the Einstein static universe. 

More can be concluded. The generalisations are symmetric with respect 
to the Einstein static universe. By symmetric one means the following. 
Instead of writing the pressure as a function of a depending on the given 
values poi Pc and A one can eliminate the cosmological constant with the 
pressure at the second centre P{a = vr), given by (|3.53|) . Let Pd and Pc2 
denote the pressures of the first and second centre, respectively. Then the 
pressure can be written as 

2PclPc2jz + {Pcl+Pc2) + iPcl-Pc2)cOSa ^ ^ 

P{a)[PcuPc2] = Po ^ . . 3.59 

2po + (Pel + Pc2) - [Pel - Pc2) cos a 



This implies 



P{^ + a)[Pci,Pc2]=P{^-a)[Pc2,Pci]. (3.60) 



Thus the pressure is symmetric to a = 7r/2 if both central pressures are 
exchanged and therefore this is the converse situation to the case where 
Ao < A < Ae. 

Another way of looking at it is to consider the pressure in terms of 
effective values, mentioned in section 1^31 Then the pressure reads 

pcffr^^ ff-l^C^l^cosa 

^ (")=^° 3-C«ffcosa' ^^-^'^ 
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where 




(3.62) 



The effective pressure at the second centre is given by — P^^{ct = vr). 
Writing P^^ for the effective pressure at the fist centre and use (|3.61|) one 
can express p^^ by the two effective central pressures. Since this leads to 
an equation quadratic in p"^ there exist two effective constant densities for 
given effective central pressures. Therefore the decomposition of the effec- 
tive energy-momentum tensor (|2.19|) in perfect fluid part and cosmological 
constant part is not unique. 

If the central pressures are equal the dependence on a vanishes and one 
is left with the Einstein static universe. 

The pressure of solutions with increasing pressure near the first centre 
is shown in figure OKI The radial coordinate is plotted in figure EHl see the 
last two cases. 
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Figure 26: 
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with Pq = 1, Pc = 1 and A = 307r 
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3.4.7 Solutions with geometric singularity 

A>As 

In this case it is assumed that A exceeds the upper hmit A5. Then (|,''{.57|) 
implies that the pressure is increasing near the centre and diverges before 
a = vr is reached. In appendix B it is shown that the divergence of the 
pressure implies divergence of the squared Riemann tensor HB.3|) . Thus these 
solutions have a geometric singularity with unphysical properties. Therefore 
they are of no further interest. 

Figure 123 shows the divergent pressure. 

Pressure 




Figure 27: Pressure function P{a) with pQ = I, Pc = 1 and A = SOvr 
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3.5 Overview of constant density solutions 
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hyperboloid 


stellar model with exterior Schwarzschild- 
anti-de Sitter solution 


A = -Snpo 


Euclidean 


stellar model with exterior Schwarzschild- 
anti-de Sitter solution 


-S-rrpo < A < 


3-sphere 


ill 11 ' 1 ^ 1 ' (HI 

stellar model with exterior Schwarzschild- 
anti-de Sitter solution 


A = 


3-sphere 


stellar model with exterior Schwarzschild 
solution 


< A < AirpQ 


3-sphere 


stellar model with exterior Schwarzschild- 
de Sitter solution 


A = Airpo 


3-sphere 


stellar model with exterior Nariai solution 


Anpo < A < Ao 


3-sphere 


decreasing pressure model with exte- 
rior Schwarzschild-de Sitter solution; the 
group orbits are decreasing at the bound- 
ary 


Ao < A < Aij 


3-sphere 


solution with two centres, pressure de- 
creasing near the first; generalisation of 
the Einstein static universe 


A = Aij 


3-sphere 


Einstein static universe 


Ae<A<As 


3-sphere 


solution with two centres, solution in- 
creasing near the first; generalisation of 
the Einstein static universe 


A> As 


3-sphere 


increasing pressure solution with geomet- 
ric singularity 
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4 Solutions with given equation of state 



This chapter analyses the system of differential equations for a given mono- 
tonic equation of state. The choice of central pressure and central density 
together with the cosmological constant uniquely determines the pressure. 
The chapter is based on jl9j . where the following analysis was done without 
cosmological term. Geometrised units where = 1/A = 1 are used. 

The existence of a global solution for cosmological constants satisfying 
A < 47r/9f, is shown, where pi^ denotes the boundary density. 

It turns out that solutions with negative cosmological constant are simi- 
lar to those with vanishing cosmological constant. For positive cosmological 
constants new properties arise. 

In section 14.11 Buchdahl variables are introduced. These new variables 
will help to proof the existence of global solutions and of stellar models. 
They are also helpful to derive analogous Buchdahl inequalities. 

Section [4.21 shows the existence of a regular solution at the centre. After 
remarks on extending the solution and a possible coordinate singularity the 
existence of global solutions is shown in section H31 

For stellar models an analogous Buchdahl inequality is derived. Solutions 
without singularities are constructed. In particular for positive cosmological 
constants this is interesting because a second object has to be put in the 
spacetime to get a singularity free solution. 

The chapter ends with some remarks on finiteness of the radius. 

4.1 Buchdahl variables 

In chapter [21 metric function v[r) was eliminated in 1)2. 7|) and ()2.8() to give 
the TOV-A equation. Alternatively the pressure can be eliminated in the 
same equations. To do this Buchdahl |2] originally introduced the following 
variables: 

C = e'^/^ 
X = . 

It seems to be natural to do exactly the same, but using the new expression 
for e~"'^'^^ which involves the cosmological constant. Alternatively it can be 
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seen as e ""^^^ written in effective values. The new Buchdalil variables are 

y2 ^ g-a(r) ^ ^ _ 2yj(^ry - -r^ (4.1) 

3 

C = (4.2) 

X = r^. (4.3) 

The derivation of the equation one looks for is unpleasant. Anyway it is a 
straight forward calculation only involving basic algebra. Nonetheless the 
derivation is shown in detail because of the interesting result: The contri- 
bution of the cosmological constant is fully absorbed by the new Buchdahl 
variable y, no other terms involving A appear in the equation. 

Start with the second field equation (|2.7|) . Putting in the new variables 
leads to 

y"^ + y^^Jx2^Jx i/^^. - 1 = (SvrP - A) x. 

With 



-2- 3, 



this gives 



8^P -\k = \y^^-2w, (4.4) 

where ^ = 2y/x-^ was used. The term ^ comes into the equation because 

i/^a; = 2^ which follows from its definition 1)4. 2|) . Differentiating with respect 
to X implies 



SvrP, = 4 2yy,,^ - y^ ^' " - 2u;,,. (4.5) 




Next the conservation equation 1)2. 8|) is rewritten. This then gives 

= -^{P + p) = -^[p + i-(3u; + 2xw,,)^ , (4.6) 

where ()2.18|) was used to eliminate p. Putting equations ()4.4|) and (|4.5j) in 
1)4. 6|1 eliminates the pressure. The equation in full details is 




2 CixC C,x \ 1 
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It is obvious that the terms with y Ci drop out. Multiplying with ^ and 
using that 

22/2/,x = -2w - 2xw^x - ^, (4.7) 
the following equation is obtained 

y'^Cxx - (w + xw^x + ^ j Cx - \'^,xC. = 0. (4.8) 



Using again 1)4. 7|) one finds that 



y^Cxx + yy,xC,x - = o, (4.9) 

or rewritten using the product rule 

fc), - = 0. (4.10) 

The last two differential equations are identical to these obtained without a 
cosmological constant and as mentioned above, A is fully absorbed by the 
new variable y. The equation is linear in C, and w. 

Derivation of (|4.10j) can also be done by rewriting the equation with 
vanishing cosmological constant in effective values for pressure and density. 
This leads to the Buchdahl variable y with Wcs- 

4.2 Existence of a unique regular solution at the centre 

It is shown that there exists a unique regular solution in a neighbourhood 
of the centre r = for each central density and given equation of state. To 
do this the following theorem is needed, see JHl for a proof. 

Theorem 1 Let V he a finite dimensional real vector space, N '.V ^ V a 
linear mapping, G : V x I ^ V a C°° mapping and g : I ^ V a smooth 
mapping, where I is an open interval in R containing zero. Consider the 
equation 

sf^+Nf = sG{s,f{s))=g{s) (4.11) 

for a function f defined on a neighbourhood of and I and taking values in 
V. Suppose that each eigenvalue of N has a positive real part. Then there 
exists an open Interval J with S J C / and a unique bounded function 
f on J\0 satisfying Moreover f extends to a C°° solution of J^. i j| ) 

on J if N,G and g depend smoothly on a parameter z and the eigenvalues 
of N are distinct then the solution also depends smoothly on z. 
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For a given equation of state equations ()2.16|) and (|2.18|) are forming 
a system of differential equations in P{r) and w{r). The aim is to apply 
theorem ^ 

Recall (|2?T8|) and (|2?T6|) with x = and dr = These substitutions 
give 

dw 3 , , s 

X— + -w = 2ttp (4.12) 
dx 2 



dp__ (dPy^ 1 {AttP + w - f ) (P + p) 



V c^P y 2 1 - 2wx - ^x 



(4.13) 



Looking for the mathematical structure of the equations one has that the 
first one is singular at the centre x = whereas the second one is regular 
there. Both equations have different properties, but using the substitution 

p = Pc + xpi, 

where pc is the value of the central density, the second equation becomes 
singular at the centre, too. The transformed system is 

dw 3 

X— + -w = 2-Kpc + 27rxpi (4.14) 
dx 2 

dpi ^ fdPy' 1 i^^P + w-^)iP + pc + xpi) 



dx \dp J 2 1 - 2wx 



X 



3' 

The system will be put in the required form to apply theorem ^ Af- 
ter substituting in p — Pf^ -\- xp\ and the corresponding pressure relation 
P = Pc + xPi{pi), only one more algebraic relation is needed. It is easily 
verified that 

/ A ( A \ f A 

( 1 — 2wx — —X I =1+1 2wx + —X 1(1 — 2wx — —x j 

Then the matrix has the form 

/ 3/2 \ 

(Pc+Pc) T • (4.16) 

\ 2dP/dp{pc) ^ J 

It is easy to obtain that N does not change if it is written in terms of effective 
values. 

Now theorem^ can be applied. Thus the system has a unique bounded 
solution in the neighbourhood of the centre, which is C°°. Therefore one 
has a unique smooth solution to (|2.16|) and (|2.18j) near the centre. 
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This disproves Collins |3| statement that for given equation of state, 
central pressure and cosmological constant the solution is not uniquely de- 
termined. 

Standard theorems for differential equations imply that the solution can 
be extended as long as the right-hand sides of (|2.1Hj) or (|4.1H|) are well 
defined. If the pressure satisfies P < oo and the denominator of (|4.13|) does 
not vanish, this means y = 1 — 2wx — (A/3)x > 0, then the right-hand sides 
are well defined. The second term involve the cosmological constant and 
therefore some new features will arise. 

Uniqueness of the solution at the centre implies the following theorem. 

Theorem 2 Suppose p{P), Pc an-d the cosmological constant A are given 
such that 

4^P, + - I = 0, (4.17) 

note that (47r/3)p(Pc) = w^c- Then the solution is the Einstein static universe 
with Ae = A. 



4.3 Extension of the solution 

Theorem 3 Assume the pressure is decreasing near the centre, this means 

Att a 
AttPc + YP(^c) - 3 > 0. (4.18) 

Then the solution is extendible and the pressure is monotonically decreasing 
if 

AttP + w-A/3>0. (4.19) 

Proof. Suppose that p = p{P), Pc and A are given such that P is decreasing 
near the centre, then w{x)^x ^ 0. Using y > because of the metric's 
signature, 1)4. 1U() implies 

(2/Cx),x < 0. (4.20) 

Rewriting ()4.4() gives 

yC,x = ^(^^7TP + w-^y (4.21) 



46 



next using the implication of H4.20() leads to 



yC,x < {yC,x) (0). 

Together with the explicit expression of (|4.21|1 this finally shows 

inP + w- ^ , , 

y > IT- 4.22 

Therefore the new Buchdahl variable y cannot vanish before the numerator 
does. Thus the right-hand sides of (|2.16|) or (|4.1.S)) are well defined and one 
can extend the solution if 4ttP + w — A/3 > 0. 

Since y > and AnP + tt; — A/3 > the sign of the right-hand side 
of (|4.13|1 is strictly negative. Therefore the density and by the equation of 
state the pressure are monotonically decreasing functions. 



4.4 A possible coordinate singularity y = 

For some cosmological constants the new Buchdahl variable y may vanish 
before the pressure does. In the constant density case the vanishing of y 
corresponded to a coordinate singularity and one could extend the solution. 
Therefore one would like to show that this is also true for a prescribed 
equation of state. 

The following theorem only indicates that the vanishing of y corresponds 
to a coordinate singularity. The square of the Riemann tensor in Buchdahl 
variables (|B.6|) is given in appendix B. 

Theorem 4 The squared Riemann tensor iB.f^) does not diverge as y — > 
if the pressure is finite. 

Proof. This is seen by first showing that y — > implies AttP + w — A/3 — > 0, 
which is the numerator of the TOV- A equation (|2.16j) and of (|4.22j) . 

Assume that p{P), Pc and A are given such that the pressure is decreasing 
near the centre. Assume further that A is given such that y vanishes before 
the pressure. Then there exists > such that 47rP<j + Ws — A/3 = 0. 
Show that Tjj = r, where f denotes the radius where the TOV-A equation 
becomes singular. This is proved indirectly: 

Assume 7^ r, i.e. ys 7^ 0. Then 47rPs + Wg — A/3 = implies = 
and equation (|4.21j) gives 

Cecils) = 0. (4.23) 
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Evaluating ()4.9() at Xg leads to 

yXxxixs) - ^C{xs)w,x{xs) = 0. (4.24) 

The pressure is assumed to be decreasing near the centre, thus w^xi-Cs) ^ 0. 
Since / (|4.24j) implies 

Cxxi^s) < 0. (4.25) 

Using 1)4. 2|) . the definition of implies that P"{rs) > 0. On the other 
hand the above assumptions on the pressure are such that P"{r) is smaller 
than zero near the centre. Therefore P"{r) < on < r < r^. And both 
inequalities together imply P" [rg) = 0, which using (|4.2|) combines to 

C,xx{xs) = 0. (4.26) 

Evaluating at Xg using ijOH]) and (jl^ finally gives 

Q{xs)w^x{xs) = 0. (4.27) 

Since C, cannot vanish w^x{xs) = 0. It was assumed that rs-,Xs > thus 
equation (|4.12j) reads 

, , 3 /47r \ 

= XsW^x[Xs) = 2 ( ~3"'^'' - I , (4.28) 

which gives Wg = (47r/3)/?s. Therefore p{r) = ps = const, on < r < r^. 
Conversely P{r) = const., which contradicts the above assumption that the 
pressure is decreasing near the centre. Thus -P'(r<j) ^ 0, therefore y has to 
vanish and = r. So the situation is similar to the constant density case. 
Therefore the expression 

Cx 1 / ^ttP + w- 4\ 

does not diverge as y — > and shows finiteness of the squared Riemann 
tensor. 

The first term of (|B.6|) gives 4/x^, the second and the third term vanish 
because ()4.29() implies y^C.x/C = 0- By the above the last term is bounded 
if the pressure is finite. 

This indicates that y = corresponds to a coordinate singularity. 

At the possible coordinate singularity the pressure has a well defined 
value given by 

P{r) = ^ - ^(O) • (4.30) 
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4.5 Existence of global solutions with A < Anpi, 



The following theorem shows the existence of global solutions. Two possi- 
bilities occur. Either the matter occupies the whole space or the pressure 
vanishes at a finite radius, in which case a vacuum solution is joined. 

Theorem 5 Let an equation of state be given such that p is defined for 
p > 0, non-negative and continuous for p > 0, C°° for p > and suppose 
that dp /dp > for p > 0. Further let the cosmological constant he given 
such that A < Airpi,. 

Then the pressure is decreasing near the centre and there exists for any 
positive value of central pressure Pc a unique inextendihle static, spherically 
symmetric solution of Einstein's field equations with cosmological constant 
with a perfect fluid source and equation of state p{P)- 

-(/A < then the matter either occupies the whole space with p tending to 
Poo CLS r tends to infinity or the matter has finite extend. In the second case 
a unique Schwarzschild-anti-de Sitter solution is joined on as an exterior 
field. 

// < A < 47rpfc the matter always has finite extend and a unique 
Schwarzschild-de Sitter solution is joined on as an exterior field. 

Proof. If the cosmological constant is given such that A < 47r/3;, then 



and the pressure is decreasing near the centre, see (|4.18|) . 

Since the pressure is decreasing near the centre the denominator of (|4.22|) 
is a positive number T>. Then the numerator of (|4.22|) can be estimated by 



in A 




y> 



V 



> 




(4.31) 



V 



> 




(4.32) 



Thus if 



A < 47r/3b 



(4.33) 
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then^ the new Buchdahl variable cannot vanish before the pressure, pf, is 
given by the equation of state because ph = p{P = 0). The coordinate Xh 
where the pressure vanishes wih be taken as the definition of the stehar 
object's radius r^,. 



If A < the matter can occupy the whole space because 1)4. 22|) implies 
positivity of y. 

Suppose P{xb) = 0. At the corresponding radius the Schwarzschild- 
anti-de Sitter solution, (|2.22)) with negative cosmological constant, is joined 
uniquely by the condition M = m(rfe). In this manner the metric is only, 
because the density at the boundary may be non-zero. The metric is at 
P{i~b) = if Gauss coordinates relative to the hypersurface P{rh) = are 
used. These are given by xi^") = Jq e^^^^/^ds. If the boundary density does 
not vanish the Ricci tensor has a discontinuity. Thus the metric is at most 
C^. Without further assumptions on the boundary density this cannot be 
improved. 

Assume now that P{x) > for all x > 0. P{x) is monotonically decreas- 
ing, therefore lima;^oo P{x) = Poo exists. This implies that the pressure 
gradient tends to zero as x ^ oo. Because of — > as x ^ oo equa- 
tion (|2.16jl does not imply that Pqo = 0, which it does if A = 0. Thus the 
equation of state only gives p — > poo = p{Poo) as x — > oo. 

< A < 47rpfe 

If < A < 4:-Kpiy then one can estimate the pressure (|4.30() at the possible 
coordinate singularity y = 0. Note that 



Hence there exists rf, such that P{rb) = 0. Since the pressure is decreasing 
and i-'(r) < it follows that r^, < f. 

Thus if the cosmological constant is positive and A < iirpi, then the pres- 
sure always vanishes at some x^. At the corresponding the Schwarzschild- 
de Sitter solution is joined uniquely by the same condition M = m(rb). The 

^The argument of |14| is applicable if the cosmological constant satisfies this condition. 



A < 




(4.34) 



which holds for all r. Therefore 




(4.35) 
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metric is at most C because the boundary density is larger than zero, this 
cannot be improved. 

4.6 Generalised Buchdahl inequality 

For stellar models an analogous Buchdahl inequality is derived. It turns out 
that the new inequality nearly coincides with the former (|3.34() . which was 
derived for constant density. 

The proof of the following theorem solves the field equations with con- 
stant density written in Buchdahl variables. This solution is compared with 
a decreasing solution and leads to an analogous Buchdahl inequality. 

The solution with constant density is denoted with a tilde above. 

Theorem 6 Let the cosmological constant be given such that A < Airpi,. 
Then for stellar models there holds 



Proof. Assume that p{P), Pc and A are given such that the pressure is 
decreasing near the centre. In section it was shown that this means that 
pressure and mean density are decreasing functions. 
Equation H4.1(J() implies 



where D is a constant of integration. 

If the density is constant C, is the only unknown function in Einstein's 
field equations. 

The function Q can be normalised to give Cfe = Ub- The constant D is 
obtained from (|4.4j) evaluated at the boundary, therefore 




(4.36) 



(4.37) 




(4.38) 



which can be used to integrate ()4.37|) . Notice that (|4.7() reads 



-2yy,x = 2w + -. 



The right-hand side of (|4.37j) can be written as 




) 



-1 
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Substitute in —2yy^x because of the first factor and divide by y. Integration 
then leads to 

~ ~ 4 / \ 

C[x) = jry[x) + const. 

The constant of integration is 

const. = C(0) + :r—\- 
2w + ^ 

So the solution of differential equation (|4.,'-{7|1 is obtained. Equation (|4.2()|1 
implies 

K,x = {yC,x)h < yC,x- 

Using that y > y one finds 

w- - ~ w - - 

C{x)>C{x) = --—^y{x) + m+ ' 



Evaluating this at the boundary and using the normalisation condition it is 
found that 

Since C(O) is positive some algebra leads to the analogous Buchdahl inequal- 
ity 

One compares solutions with decreasing mean density and solution with con- 
stant density. The compared constant density corresponds to the boundary 
mean density. Therefore (|4.39j) reads 

which nearly equals the analogous Buchdahl inequality 1)3.34(1 if the expres- 
sion for is used. The only difference is that the equation contains the 
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boundary mean density Wh rather than the boundary density pf,. Exphcitly 
written out leads to 



l-2u;fer2--r2> - (4.41) 



A 1 A 

—rr > 

3^-3 9wb 

and holds for all monotonically decreasing densities and proves theorem |SJ 



Equation H4.34() implies that the right-hand sides of the last two equa- 
tions (|4.4H) and (|4.4flj) are positive. This improves (|4.22j) . which only gave 
y > 0. Solving inequality ()4.41|) for r;, gives an analogue of 1)3.35(1 

9 9 V 3i«i, y , , 

Tb < — —■ 4.42 

As before the only difference is that the boundary mean density is used 
instead of the boundary density itself. Multiply l|4.42() by Wb and use that 
Wb = m{rb)/r^ gives 



rb 9 \ 3'm{rb) 
Taking all terms on one side and multiplying by 9m{rb)rb leads to 



9m{rb)'^ - Am{rb)rb + < 0. 



Rewriting this with 



2 \ 2 ^ 

9m{rbf - Am{rb)rb = ( 3m(rb) - -n] - -r^, 



finally gives 



3m{n) - '^n) - rl - ] < 0. (4.43) 



Another way of writing the above (j4.43j) is 



3m(rb) - ^rfe - n^J ^ - ^rfj ^3m(rb) - + n^J ^ - ^rf j < 0, 



(4.44) 
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where a? — l)^ = {a — b){a + h) was used. Putting A = shows that the second 
term has to be positive because of positivity of the mass. Then the first 
term is negative and imphes the Buchdahl inequahty (|3.38() . For negative 
cosmological constants the square root terms are always weh defined. For 
positive cosmological constants they are well defined if 

n < (4.45) 

Then both expressions of the product (|4.44j) can be combined to 



With A = equation (|4.46|) directly implies Buchdahl's inequality (|H.H9j) . It 
holds for arbitrary static fluid balls in which the density does not increase 
outwards. 



4.7 Solutions without singularities 

The last two sections showed the existence of stellar models for cosmological 
constants A < ^i^pb- Equation (|4.22|) implies that the boundary of the 
stellar object has a lower bound given by the black-hole event horizon and 
an upper bound given by the cosmological event horizon. The upper bound 
only occurs if the cosmological constant is positive, as already said. 

Stellar models with A < have a lower bound given by the black-hole 
event horizon. At the boundary the Schwarzschild-anti-de Sitter solution 
is joined on as an exterior field. Figure [29l shows that this joined exterior 
field has no singularities. Therefore stellar models with A < have no 
singularities. Solutions with cosmological constant satisfying A < are 
globally static. 

For positive cosmological constants the situation is different. But one 
can construct solutions without singularities. 

At the boundary r = the Schwarzschild-de Sitter solution is joined 

by the usual procedure introducing Gauss coordinates. This leads to 
Penrose-Carter figure QUI This spacetime has an infinite sequence of singu- 
larities r = and space-like infinities r = oo. 

The surface r = r;, can also be found in the vacuum region where the 
time-like Killing vector is past directed. This means that a second stellar 
object is put in the spacetime. Therefore Penrose-Carter figure mi gives a 
general picture and is not restricted to constant density solutions. 
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Note that solutions with positive cosmological constant are not globally 
static because of the dynamical parts of the exterior Schwarzschild-de Sitter 
solution. 

Solutions without singularities and without horizons were recently de- 
scribed by They considered an interior de Sitter region and an exterior 
Schwarzschild solution separated by a small shell of matter with equation of 
state p{P) = P. The shell replaces both the de Sitter and the Schwarzschild 
horizon. The new solution has no singularities. 



4.8 Remarks on finiteness of the radius 

So far it has been shown that given an equation of state, a central pressure 
and a cosmological constant there exist a unique model with finite or infinite 
extend. This depends on the given equation of state and on the cosmological 
constant. This section gives criteria to distinguish these two cases, see jl9| . 



A > 

If A > solutions are always finite. Thus given an equation of state p = p{P) 
and cosmological constant such that 

A < 4ttp{P = 0) = Airpb, (4.47) 

then there always exist a radius rb where the pressure vanishes. 



A < 

If A < either the pressure vanishes for some finite radius or the density is 
always positive and tends to poo as r tends to infinity. 

First a necessary condition for finiteness of the radius is derived. Recall 
1)2. 16(1 written in the variable x 

1 UttP + w-^)(P + p) 
Px = — ^ —■ 4.48 

Assume that the stellar model is bounded by Xf, or the corresponding r^- 
On the closed interval [0, Xb] the function y has a positive minimum K. The 
maximal contribution due to the first factor is given by its value at the 
centre. Therefore 



>-L{P + p), 
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where the constant L is given by 



L=^. i^^Pc + wc-^]. (4.49) 



2A' V 3 
Compare the above inequahty with 

which can be integrated to give 

dP 



P + p 



Thus we obtain a bound for the integral. Hence if the star has finite bound- 
ary then 

f^" dP 

/ _<oo. (4.50) 

This gives a necessary condition for the equation of state to give fluid balls 
of finite extend. Then equation 1)2. 17(1 implies that z^(r) is finite. 

Now a sufficient condition is derived. Since the density distribution is 
decreasing m > ^pr^, w > for all r > 0. 

Thus (|4.48j) gives 

1 27r n , , 

P,x<-^wp<-—p^, (4.51) 



where one needed that 



_ I > 0, 

which is always achieved as long as A < 0. Then comparing with a P* again 
and using that there has to exist a point Xb where the pressure vanishes one 
may integrate ()4.51|) . Therefore if 

^'^ dP 

^ < oo, (4.52) 

P 

the stellar object has finite radius. Both criteria (|4.5nj) . (|4.52|) only involve 
the low pressure behaviour of the equation of state. 
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Appendices 

A Einstein tensor and energy- momentum conserva- 
tion 

In the following a static, spherically symmetric metric is considered. In the 
used notation it may be written as 

ds"^ = -\e^''('^Ut'^ + e'^'^'-Ur^ + r\de^ + sin^^ dc/)^). (A.l) 
A 

The non vanishing ChristofFel symbols are: 

rg, = -re-'^W r;^ = -rsin2^e-'^W 

Tqj. = - = — cos 9 sin 9 

^tr = l rJ, = COt0. 

The trace terms are: 

rf. = r,-^ = ^+l(A,.'(r)+a'(r)) 

r^^ = cot0 r5. = o. 

The Ricci tensor is defined by 

Its components are: 

= e^-^W-^W ( -u"(r) + -Xu'(rf + -i/'(r) - -a'(ry(f 
\2 4 r 4 ' 

i?,, = -lxi^"{r) - h^i^'irf + ia'(r)Ai/'(r) + -a\r) 
2 4 4 r 

Rfjxl) = sin^ ^'-Ree 
= if /_i / I/. 
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The Ricci scalar is 

The Einstein tensor is defined by 

Only its four diagonal components remain. They are given by: 

Gtt = Ae^'^'^ ^[^- ^e-"^'^) (A.4) 

C = \(\^T\v\r)-e<'^\ (A.5) 



GeQ = r2e-'^W^(Az."(r) - -a'(r) + -Xi^'{r) + ^A2i/'(r)2 - ia'(r)Ai/'(r)) 
2 r r 2 2 

(A.6) 

G^^ = sin2 ^Gee. (A. 7) 

The Einstein tensor has vanishing covariant derivative and therefore implies 
energy-momentum conservation, 

V^Tf"" = 0. (A.8) 

For a perfect fluid in a static, spherically symmetric spacetime the energy- 
momentum tensor has the form 

V = + XP{r))U^U, + P{r)g^u, (A.9) 

where 4- velocity Ut = {—1/X)e^'^^'^^^'^ is given for a fluid at rest. Conservation 
of this quantity gives four equations of motion. Because of symmetries, only 
the radial component /i = r does not equal zero. By using the Christoffel 
symbols derived above, it is found that 

= {p{r) + AP(r)) ^e-'^^'^) + P'(r)e-'^W, 

which gives 

2P'{r) = -{XP{r) + p{r)y{r). (A.IO) 
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B Geometric invariants 

Constant density solutions 

Metric 1)3 .21) of the constant density solutions for A > —Sirpo can be written 
ds^ = -( -p^^^ ' dt^ + i(da2 + sin^a dn^). (B.l) 
Integration of gives 

P{c^) = Po ^ \ ' . B.2 

3 — C cos Q 

From this one can find the square of the Riemann tensor. It reads 

2 

V 3 '^^ ' 3 



81 - C72 cos2 a (9 - 6C7 cos a - iC^ cos^ q) 
(3 + C cos Oif' (3 — C cos a)^ 



(B.3) 



and gives two important implications. 

First it shows that the radial coordinate r behaves badly as r — > r, which 
corresponds to a ^ 7r/2. Thus 



R,,uXnR^''^\a = n/2) = 12 (^^po + |) , (B.4) 

and one clearly has a coordinate singularity because the square of Riemann 
tensor is finite. It was already pointed out that this is only a coordinate 
singularity when the geometry of the interior solutions was discussed. 

Secondly, the square of the Riemann tensor diverges as the pressure 
diverges. Compare the denominators of (|B.2|) and HB.3|) . Thus there is a 
geometric singularity if the pressure is divergent. 

Buchdahl variables 

In Buchdahl variables 



= e 



2 _ „-a{r) ^ ^ _ 2w{ry - -r^, 



vl2 2 

e ' X = r , 
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the static spherically symmetric metric can be written 



-Qixfdt^ + 



+ x{d9 + sin"^ 



4xy(x)2 

In these coordinates the square of the Riemann tensor becomes 



(B.5) 



+ 2 



2(1/^), 



+ 



(B.6) 



The square of the Weyl tensor is given by 



Conformal flatness of constant density solutions 



(B.7) 



In chapter 0] Einstein's field equation with cosmological constant for a perfect 
fluid were rewritten in terms of Buchdahl variables. This gave equation 
(|4.10j) . stated again 



1 •w.xC 

2 y 



0. 



(B.8) 



If the density is assumed to be constant then w^x = 0. Therefore equation 
(|B.8|) implies that the square of the Weyl tensor (|B.7|) vanishes. In spherical 
symmetry this then implies that the Weyl tensor vanishes. Thus constant 
density solutions are conformally flat. Furthermore it follows that this is 
not true if w^x 7^ 0. 

Invariants for other metrics 

To complete this appendix the invariants of the Schwarzschild-de Sitter 

1)2. 22() and Nariai (|3.48|) solution are given. For Schwarzschild-de Sitter 
they are given by 



Rf^.x.R^"^^ = 48— + -A2 



r" 3 

M2 



(B.9) 
(B.IO) 
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The Nariai metric gives 

^^.Ak^^"^" = SA^ (B.ll) 
C^.AkC^'^^" = ^A2. (B.12) 
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C Penrose-Carter diagrams 




Figure 28: Penrose-Carter diagram for Schwarzschild-de Sitter space. The 
Killing vector d/dt is time-like and future-directed in regions / and time-like 
and past-directed in regions IV . In the others regions it is space-like. The 
surfaces r = r_|_ and r = r_|__(_ are black- hole and cosmological event horizons, 
respectively. J^"*" and J"" are the space-like infinities. 




Figure 29: Penrose-Carter diagram for Schwarzschild-anti-de Sitter space. 
The surface r = r+ is the black-hole event horizons. J^"*" and are the 
time-like infinities. 
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Figure 30: Penrose-Carter diagram with one stellar object having a radius 
R which lies between the two horizons. The group orbits are increasing 
at the boundary. The shaded region is the matter solution with regular 
centre. There is still an infinite sequence of singularities r = and space- 
like infinities r = oo. 



r = oo, J^"*" 




r = oo, 



Figure 31: Penrose-Carter diagram with two stellar objects having radii R 
which lie between the two horizons. Since the group orbits are increasing 
up to R the vacuum part contains the cosmological event horizon r++ . This 
solution with two objects has no singularities. Because of regions He and 
IIIc this spacetime is not globally static. 
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Figure 32: Penrose-Carter diagram with two stellar objects having radii 
R which lie between the two horizons. The group orbits of the interior 
solutions are decreasing where the vacuum solution is joined. Thus the 
r = singularity of the vacuum part is present. 
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Deutsche Zusammenfassung 



Die vorliegende Arbeit untersucht statische, kugelsymmetrische Losungen 
der Einsteinschen Feldgleichungen mit kosmologischer Konstante. Fiir die 
Materie wird eine ideale Fliissigkeit mit vorgegebener monotoner Zustands- 
gleichung angenommen. Es werden neue Losungstypen beschrieben. 

Bemerkenswert ist, daB die kosmologische Konstante innerhalb dieser 
Fragestellung lange vernachlassigt wurde. 

Das erste Kapitel ist eine Einfiihrung. Sie ordnet die Thematik innerhalb 
der Gravitationstheorie ein und geht auf Teile vorhandener Literatur ein. 

Obwohl schon Weyl |2Z| diese Losungen untersuchte und wesentliche 
Ergebnisse darstellte, wurden sie kaum weiter beachtet. Erst wieder Kriele |13j . 
Stuchlik und Winter untersuchten solche Losungen. 

Im zweiten Teil werden die relevanten Gleichungen hergeleitet und deren 
mathematische Struktur beschrieben. 

Wird fiir die Materie eine ideale Fliissigkeit angenommen, so erhalt 
man eine Verallgemeinerung der Tolman-Oppenheimer-Volkoff Gleichung 
mit kosmologischer Konstante ()2.16|) . kurz TOV-A Gleichung. Sei nun eine 
Zustandsgleichung p = p{P) vorgegeben. Dann bilden die beiden Gleichun- 
gen ()2.16|) und (|2.18() ein System von gewohnlichen DifFerentialgleichungen 
in den Funktionen P(r) und w{r). 

Fiir das Vakuum lassen sich die Feldgleichungen einfach losen. Man 
erhalt die bekannte Schwarzschild-de Sitter Losung (|2.22j) fiir positive kos- 
mologische Konstante bzw. die Schwarzschild-anti-de Sitter Losung fiir 
das negative Vorzeichen. Dies ist die einzige statische, kugelsymmetrische 
Losung der Vakuumfeldgleichungen mit nichtkonstanten Gruppenbahnen. 

Penrose- Carter Diagramm |2H1 zeigt die Schwarzschild-de Sitter Raum- 
Zeit mit dem Ereignishorizont des Schwarzen Loches und dem kosmologis- 
chen Ereignishorizont. 

Das Penrose-Carter Diagramm 1291 zeigt die Schwarzschild-anti-de Sitter 
Raum-Zeit. In ihr gibt es nur den Ereignishorizont des Schwarzen Loches. 

Mit kosmologischer Konstante gibt es eine weitere statische, kugelsym- 
metrische Losung der Feldgleichungen fiir das Vakuum, die Nariai Losung 
|16| I17j. Das Volumen der Gruppenbahnen ist hierbei konstant. Finer der 
neuen Losungstypen erfordert diese Losung fiir den Vakuumteil. 

Nach dieser Beschreibung der Gleichungen wird kurz auf deren Newton- 
schen Grenzwert eingegangen. Wie zu erwarten, ergibt sich die Newtonsche 
Theorie mit einem zusatzlichen A- Term (|2.23|) . Das Newtonsche Analogon 
der TOV-A Gleichung ist die Euler Gleichung der Hydrodynamik im hydro- 
dynamischen Gleichgewicht, wieder mit korrigierten Potential, 1)2. 25|) . 
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Der dritte Abschnitt befaBt sich mit den speziellen konstante Dichte 
Losungen. Die TOV-A Gleichung wird fiir den einfachen Fall einer Zus- 
tandsgleichung p = p{P) ^ p = po = const., namlich homogener Dichte, 
integriert. Der zusatzlich vorhandene A- Term eroffnet neue Losungen. Die 
t = const. Hyperflachen sind bei verschwindender kosmologischer Konstante 
3-Spharen. A kann nun auch so gewahlt werden, daB diese Hyperflachen 
euklidisch oder hyperbolisch werden. Die Losungen mit konstanter Dichte 
sind konform flach, unabhangig von der kosmologischen Konstanten. Mit 
den im vierten Kapitel eingefiihrten neuen Buchdahlvariablen kann leicht 
gezeigt werden, daB dies direkt aus den Feldgleichungen folgt. 

Konstante Dichte Losungen wurden nur bis zu einer oberen Grenze 
A < 47r/3o untersucht [21]. Bis zu dieser Grenze beschreiben die Losun- 
gen ausschlieBlich Modelle von Sternen. Es existiert immer ein Radius, an 
dem der Druck verschwindet und an dem eine Vakuumlosung angeschlossen 
wird. Fiir diese Modelle laBt sich eine analoge Buchdahl Ungleichung her- 
leiten. Sie besagt, dafi der Radius stets zwischen dem Ereignishorizonts des 
entsprechenden Schwarzen Loches und dem kosmologischen Ereignishorizont 
liegt. Der kosmologische Ereignishorizont existiert nur im Falle positiver 
kosmologischer Konstante. Daher haben Modelle mit A < im allgemeinen 
keine obere Schranke fiir den Radius. 

Diese Losungen konnen frei von Singularitaten konstruiert werden. Falls 
A < ist dies offensichtlich, da die Aufienraumlosungen frei von Singu- 
laritaten ist. Fiir A > muB ein zweites Objekt in die Raum-Zeit gesetzt 
werden, damit diese Losungen frei von Singularitaten werden, siehe Penrose- 
Carter Diagramm l31l 

Wird eine groBere kosmologische Konstante vorgegeben, so wird die Ko- 
ordinatensingularitat erreicht, noch bevor der Druck verschwindet. An ihr 
divergiert der Druckgradient, der Druck jedoch bleibt endlich. Da es sich 
lediglich um eine Koordinatensingularitat handelt, konnen entsprechende 
Losungen fortgesetzt werden. In diesem dritten Kapitel werden die ver- 
schiedenen Losungen durch die kosmologische Konstante parametrisiert. 

Entspricht A der Grenze, so verschwindet der Druck bei der Koordi- 
natensingularitat. Nur in diesem Fall erfordert die Innenraumlosung, dafi 
die Nariai Metrik angeschlossen wird. In den anderen Fallen ist es immer die 
Schwarzschild-de Sitter Losung, die den Vakuumteil beschreibt. Wird der 
kosmologische Term nun grofier gewahlt, dann verschwindet der Druck nach 
der Koordinatensingularitat. Dort ist das Volumen der Gruppenbahnen fall- 
end und man hat den Teil der Schwarzschild-de Sitter Metrik anzuschliefien, 
der die r = Singularitat enthalt, siehe Penrose-Carter Diagramm l32l Eine 
weitere VergroBerung fiihrt dazu, dafi der Druck nicht mehr verschwinden 
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kann, er ist strikt positiv. Es ergeben sich Losungen mit zwei Zentren, die in 
Fall konstanter Dichte beide regular sind. Dies sind Verallgemeinerungen des 
statischen Einstein Universums. Sie haben homogene Dichte aber inhomoge- 
nen Druck. Der im ersten Zentrum vorgegebene Zentraldruck fallt monoton 
bis zum zweiten Zentrum. Bei weiterer VergroBerung ergibt sich der Ein- 
stein Kosmos selbst, um dann deren Verallgemeinerungen mit ansteigendem 
Druck bis zum zweiten Zentrum zu beschreiben. In Abschnitt 1114.61 wurde 
gezeigt, dafi sich die Verallgemeinerungen des statischen Einstein Univer- 
sums mit ansteigendem Druck im ersten Zentrum und diese mit fallendem 
Druck im ersten Zentrum ineinander iiberfiihren lassen. 

Als letztes, bei noch grofierer kosmologischer Konstante ist der Druck 
monoton ansteigend und divergiert. Dort hat die Raum-Zeit eine geometrischen 
Singularitat. Diese Losungen sind unphysikalisch und werden nicht welter 
diskutiert. Am Ende dieses Kapitels werden alle konstante Dichte Losungen 
in einer Ubersicht dargestellt. 

Das vierte Kapitel beschaftigt sich zunachst mit der Existenz von glob- 
alen Losungen. Mit neuen Buchdahlvariablen kann diese bis zu einer oberen 
Grenze der kosmologischen Konstante (|4.H,''{|) bewiesen werden. Die entsprechen- 
den Losungen werden beschrieben und mit denen konstanter Dichte ver- 
glichen. Die obere Grenze, bis zu der die Existenz einer globalen Losung 
gezeigt werden kann, schlieBt aber die neuen Losungstypen aus. 

Fiir endliche Losungen wird die verallgemeinerte Buchdahlungleichung 
1)4. 41() hergeleitet. Sie gilt fiir alle statischen, kugelsymmetrischen Fliissigkeit- 
skugeln deren Dichte nach aufien nicht ansteigt. Damit wird gezeigt, dafi 
Losungen ohne Singularitaten nicht nur im Fall konstanter Dichte konstru- 
iert werden konnen, sondern dafi diese Losungen fiir beliebige Zustandsgle- 
ichungen diese Eigenschaft besitzen. 

Den AbschluB des Kapitels bilden Bemerkungen iiber die Existenz eines 
endlichen Radius fiir eine vorgegebene Zustandsgleichung. Es werden eine 
notwendige und eine hinreichende Bedingung hergeleitet. 

Im Anhang dieser Arbeit werden Ricci und Einstein Tensor sowie einige 
geometrische Invarianten berechnet. Dariiber hinaus sind dort die Penrose- 
Carter Diagramme dargestellt. 
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